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ABSTRACT 

We investigate the role of thermal fluctuations and of the finite number of monomers in 
small clusters of carbon atoms on the nucleation rate of carbonaceous grains. Thermal 
fluctuations are due to the quantized nature of the energy exchanges between the 
clusters, the gas, and the radiation field. Nanoscale effects modify the spontaneous 
detachment of monomers due to the finite amount of internal energy contained in 
small clusters. We find that both corrections have a big impact on the stability of the 
clusters and on the rate of nucleation. We implement our model within a Monte Carlo 
code to derive the new stability conditions for clusters as well as nucleation rates. 
Due to computing limitations, we can explore the consequences of this approach only 
at high temperatures, at which particle interactions are not much less frequent than 
photon interactions. We found that the combined effect of the detachment correction 
and the temperature fluctuations produces faster nucleation. We also found that the 
nucleation rate depends on the composition of the gas and not only on the partial 
pressure of the compound that condensates into grains. This is a unique result of this 
model that can be used to prove or disprove it. 
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1 INTRODUCTION 

Dust is one of the main constituents of the interstellar 
medium (ISM), making up to 1 per cent of the baryonic mass 
of galaxies that is not in stars (Mathis, Rumpl & Nordsieck 
1977). The average composition of interstellar dust is made 
of Polycyclic Aromatic Hydrocarbons (PAHs), carbonaceous 
grains, and silicates (Weingartner & Draine 2001), each of 
which contains particles ranging from several angstroms to 
a few microns in size. Dust plays many important roles in 
astrophysics: it provides attenuation of radiation (Cardelli, 
Clayton & Mathis 1989), induces polarization in the optical 
and UV frequencies (Serkowski, Mathewson & Ford 1975), 
provides a safe environment for molecular chemistry (Brown 
1990), cools protostellar clouds allowing for the formation 
of population II. 5 stars (Schneider et al. 2006), and, finally, 
provides the building blocks for the formation of terrestrial 
planets (Lissauer 1993). 

Despite such fundamental roles played by dust, the 
physics of dust nucleation is still highly debated, and the 
sites, rates, and composition of the dust produced in the dif- 
ferent environments are poorly known. Historically, dust nu- 
cleation has been treated as an out of equilibrium perturba- 
tion of a thermodynamical process (Becker & Doring 1935; 
Feder et al. 1966; Kashchiev 2000). Under this approach, the 
dust /vapor mix is assumed to be in thermodynamical equi- 



librium and able to exchange energy in a continuous way. 
A critical cluster size can be identified as the cluster with 
the highest Gibbs free energy and the nucleation rate is com- 
puted as the flow of particles through the critical cluster size 
(Becker & Doring 1935; Feder et al. 1996; Kashchiev 2000). 
Even though this "classical" theory can give an explanation 
to the qualitative features of nucleation (the presence of a 
threshold below which nucleation is forbidden, the existence 
of supersaturation, and the steep increase of the nucleation 
rate with both temperature and saturation), the quantita- 
tive comparison between prediction and rates is still prob- 
lematic (Wolk & Strey 2001). To make things worse, any 
attempt to improve the theory has resulted in a deteriora- 
tion of the quantitative comparison. 

The adoption of the thermodynamical approach for as- 
trophysical dust nucleation was criticized already by Donn 
& Nuth (1985) who pointed out problems with the time 
scale for obtaining thermodynamic equilibrium as well as 
discrepancies between the theory and experiments. An al- 
ternative to the thermodynamical treatment is the kinetic 
treatment of nucleation (Becker & Doring 1935, Kashchiev 
2000). The kinetic theory derives nucleation rates by an- 
alyzing the accretion and ejection rates of monomers from 
sub-critical clusters. The most fundamental difference is that 
in the kinetic theory the nucleation rate depends on the 
properties of grains of all sizes (and especially of the unsta- 
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ble ones) while in the thermodynamic theory the nucleation 
rate depends only on the properties of the critical cluster. 
Another advantage of the kinetic theory is that of allow- 
ing the consideration of radiation as a player in the nu- 
cleation game, as both an energy sink (Lazzati 2008) and 
as an energy source (this paper). Astrophysical dust nucle- 
ation is more complex than lab nucleation (e.g. of water, 
Wolk & Strey 2001) for at least two reasons. First, the con- 
ditions of nucleation in astrophysics are very different from 
lab conditions and the verification of any theory that pre- 
dicts nucleation rates of a particular dust compound can be 
performed only very indirectly. Second, some of the com- 
ponents of interstellar dust, such as silicates, do not have 
a vapor phase and the nucleation process must include the 
chemical bonding of the atoms into the basic monomer of 
the grain. We call this process chemical nucleation to dis- 
tinguish it from homogeneous nucleation (the nucleation of 
a vapor on its own liquid or solid phase) and from hetero- 
geneous nucleation (the nucleation of a compound on a dif- 
ferent substrate). This difficulty of chemical nucleation has 
been approached in two ways. Kozasa et al. (1989, 1991, see 
also Nozawa et al. 2003, 2008, 2010; Todini & Ferrara 2001; 
Schneider et al. 2004; Bianchi & Schneider 2007) consider 
the nucleation process as driven by the atom that accretes 
with the lowest frequency, assuming all other compounds to 
be already present on the grain surface. It is unclear, how- 
ever, how the very first molecules are formed and bound 
together in a sub-critical cluster. Cherchneff et al. (2000) 
and Cherchneff & Dwek (2009, 2010) consider instead the 
chemistry of precursor molecules to bridge the gap between 
atoms in the gas phase and sub-critical clusters onto which 
the monomers can accrete. 

In this paper we address another limitation of the ther- 
modynamical treatment by considering the effect of discrete- 
ness in the energy exchange between the clusters, the gas, 
and the radiation field. We focus on the homogeneous nu- 
cleation of carbonaceous grains to avoid the complication 
of chemical nucleation, and approximate cluster surfaces as 
spherical to simplify calculations. As a consequence of the 
quantized nature of cluster-radiation and cluster-atom inter- 
actions, the cluster vibrational temperature can grow well 
above the equilibrium temperature and/or decrease below 
it. We develop a Monte Carlo code to simulate the thermal 
history and stability of small clusters and grains and derive 
nucleation rates. Due to technical computing limitations we 
can, at present, address only the high temperature regime, 
since at low temperatures the interactions are overwhelm- 
ingly dominated by exchanges of radiation. 

This paper is organized as follows: in Section 2 we de- 
scribe the theory of kinetic nucleation and obtain equations 
for the temperature fiuctuations and ejection rates in small 
clusters, in Section 3 we describe the results of the Monte 
Carlo implementation, and in Section 4 we discuss our re- 
sults. 



of monomers (in our case Carbon atoms) from a cluster or 
grain Q 

The nucleation rate, i.e., the number of new stable 
grains formed per unit time and volume, is given by: 

-I -1 



J = nxfi 




(1) 



where nx is the density of monomers in the gas phase and 
fi is the rate of the reaction 



(2) 



where X is a monomer and Xi is a cluster that contains i 
monomers. The rate of the reaction 



Xi Xi-i + X 



(3) 



is instead indicated as gi. 

In an equilibrium situation, where the cluster and the 
gas are at the same temperature, the attachment rate is 
given by: 
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where fcs < 1 is the sticking coefficient, nx is the number 
density of the monomer in the gas/vapor phase, and px 
is the density of the solid phase. The detachment rate is 
instead given by: 



Qi = ksnx,oq{T){36TT)^^^ 
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where nx,Gq(T) is the density of the gas phase at saturation 
for the temperature T and a is the surface tension or surface 
energy of the condensed phase. The extra exponential term 
takes into account the increased detachment rate from small 
cluster due to the effect of surface tension (or surface energy) 
that decreases the effective binding energy of monomers. 

While the attachment fi only depends on the gas prop- 
erties (density and temperature), the detachment frequency 
depends on the cluster properties and can change substan- 
tially if the assumption of thermal equilibrium is relaxed. 
Following Lazzati (2008), we assume that if the grain tem- 
perature Tgrain is difi^ereut from the gas temperature Tgas, 
the detachment frequency is modified as 



Qi = fcsnx,oq(rgrain)(36-7r)^''^ ( 
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(6) 



i.e., we assume that the detachment rate depends only on 
the grain properties and not on its surroundings. Later, we 
will describe an additional modification of the gi equation to 
take into account the reduction in detachment rate from very 



2 NON-EQUILIBRIUM NUCLEATION 

The kinetic approach to nucleation (Becker & Doring 1935; 
Kashchiev 2000) allows for the computation of nucleation 
rates in any non-equilibrium condition. The basic require- 
ment is the knowledge of the accretion and ejection rates 



^ In this paper we will define a cluster as an unstable group of 
monomers, with a size smaller than the critical radius. Similarly, 
we will define a grain as a stable group of monomers, larger than 
the critical radius. However, when the stability is not the focus 
of a sentence, cluster and grain are used as synonyms. 



small and/or very cold clusters (Guhathakurta & Draine 
1989). Once Eq.[6]is established, the challenge is to evaluate 
the cluster temperature to be used. Lazzati (2008) assumed 
that the input and output of energy can be considered con- 
tinuum and found the equilibrium solution, showing that 
in most cases the cluster stabilizes at a temperature that 
is smaller than that of the gas, due to the increased radi- 
ation losses. In this paper, we improve on that assumption 
by considering the quantized nature of the energy exchanges 
between the cluster, the gas, and the radiation field. We con- 
sider five interactions responsible for altering the thermal 
state of the cluster: 

• attachment of a monomer to the cluster 

• detachment of a monomer from the cluster 

• collision between the cluster and an atom/molecule of 
an inert carrier gas 

• absorption of a photon 

• emission of a photon 

Due to the statistical nature of the occurrence of all the 
above processes, we include them in a Monte Carlo code. At 
each iteration, the code stochastically selects the event that 
would occur first. For that iteration, energy conservation 
laws are applied to determine the thermal state of the cluster 
after the selected event occurs. The temperature change is 
computed through the change of internal energy, A (7. The 
cluster size is also modified if an attachment or detachment 
event occurs. 



2.1 Monomer Attachment 

The rate of attachment of monomers is given by Eq. |31 The 
change of internal energy for a monomer attachment has 
two components: the kinetic energy gained from the incom- 
ing monomer and the binding energy that is released at the 
attachment: 

AUau = KE{Tg^,) + Esindii -^i + l) (7) 

where the kinetic energy KE{T^!is) is randomly drawn from 
a Maxwell-Boltzmann distribution at the gas temperature 
Tgaa, and i?Bind(i) = iEBind,bik " oS{i) is the total bind- 
ing energy released in the assembly of i monomers into a 
cluster with surface S(i) and binding energy EBind.uk for 
a monomer in the bulk. Thus, the binding energy for an 
attachment is computed as: 



-EBind(i -^i + l) = 

[{i + l)-BBind,bik - o-S{i + 1)] - [j£Biiid,bik - o-S{i)] = 
i?Bind,bik — a[S(i + 1) — S{i)] = 




where we have used the relation S{i) = (SGvr)^/'^ ( 'px' )^^'^' 
2.2 Monomer Detachment 

Similarly, the change of internal energy for a monomer de- 
tachment has two components: the kinetic energy lost from 
the outgoing monomer and the binding energy required to 
separate the monomer from the cluster: 
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Symbol Meaning 



b Minimum number of quanta to eject a monomer 

Cabs Correction term for Photon Absorption 

Comis Correction term for Photon Emission 

Cq Correction term for Monomer Detachment 

df Number of degrees of freedom 

^Bind Binding energy 
i?Bind,blk Binding energy of a bulk monomer 

Attachment rate of monomers onto a cluster of size i 

gi Detachment rate of monomers from a cluster of size i 

hi Collision rate of inert gas atoms with a cluster of size i 

H/C Hydrogen/Carbon number density ratio 

iQ Critical cluster size 

J Nucleation rate 

KE Kinetic energy of a monomer or inert gas particle 

ks Sticking coefficient 

mcG Mass of an inert gas particle 

mx Mass of the monomer 

ncG Number density of carrier gas 

nx Number density of monomers in the gas phase 

"x,cq Number density of monomers in the gas phase 

at vapor pressure 

P{y, i) Modified Blackbody Distribution 

q Total number of quanta in a cluster 

S Gas saturation 

S{i) Surface area of cluster containing i monomers 

Tgas Temperature of the gas 

Tgrain Temperature of the grain 

^rad Temperature of the radiation field 

U Internal energy 

X Monomer 

Xi Cluster containing i monomers 

cii Photon Absorption rate of a cluster of size i 

^Uxtt Change of internal energy for monomer attachment 

£{v, i) Emissivity, frequency and size dependent 

7i Photon Emission rate of a cluster of size i 

px Density of the condensed phase 

(T Surface energy (tension) of the condensed phase 

crsB Stefan-Boltzmann constant 

ujo Frequency of harmonic oscillators in the Einstein model 
Table 1. Glossary of used symbols 

AUoot = -A'S(Tg,ain) + Esindii ^ i - 1) (9) 



where the kinetic energy i('ii'(Tgiain) is randomly drawn from 
a Maxwell-Boltzmann distribution at the grain temperature 
Tgrain. So, the biudlug energy is: 

-EBind(i i - 1) = 

[{i - l)-E'Bind,blk - (^S{i - 1)] - [i-BBind.blk ~ (^S{i)] = 

--EBind.bik - o- [S{i - 1) - S{i)] = 

-^Bind.bik - a(36^)^/^(^^^ - 1)^/^ - t'/'] (10) 

again for 5(i) = (36^)^/«(^)'/'. 

The rate of detachment of monomers from big grains is 
be given by Eq. [5] However, as previously mentioned, gi re- 
quires a modification to correctly describe detachment from 
very small and/or very cold clusters. As Guhathakurta & 
Draine (1989) realized, clusters must have at least enough 
internal energy to supply a single particle within the grain 
with at least the binding energy to eject that particle from 
the cluster. However, even if the internal energy of the clus- 
ter is above this threshold, we cannot be certain that this 
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energy will break the bond of a single monomer. In order to 
determine the probability that this could happen, wc nmst 
introduce quantization of energy, modelling clusters as Ein- 
stein solids, such that Xi has df = 3i — 6 vibrational degrees 
of freedom. Assuming that these degrees of freedom arc har- 
monic oscillators with common frequency uio = 0.75kTDebye 
and average energy Tiujo, and given the internal energy U of 
the cluster, we can determine the total number of quanta 
in the cluster, q = U /TiuJo, as well as the number of quanta 
required to eject a monomer, b = E-Qind/T'i^n- Using a com- 
binatorial approach with the preceding values, we calculate 
the probability that one degree of freedom has at least b 
quanta. This (sec Guhathakurta & Draine 1989) reduces the 
monomer detachment rate by: 



q\{q-b + df - 1)! 
{q + df-mq-by. 



(11) 



where ^"7" + 1^ is the normalization factor for df 
Thus, the new detachment rate is simply: 

9i = CgQi 



(12) 



If we define y = q/df as the average number of quanta per 
degree of freedom, Cq is nearly unity when y ^ b while 
Cg — >■ as ?/ — ^ or g — >■ 6; consequently, for q < b, g* = Q 
and monomers cannot be ejected from the cluster. 



2.3 Inert Carrier Gas Collision 

Even though we focus on homogeneous nucleation, we also 
consider collisions with inert particles, such as hydrogen 
atoms, that are inevitably mixed in the surrounding gas in 
any astrophysical setting. We model this interaction as an 
event in which the particle exchanges energy with the clus- 
ter and rebounds with a new velocity distribution. The rate 
of inert collisions is as follows: 



hi = fesncG(367r)^^^ ( 
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(13) 



Although similar to fi, this rate is dependent on the number 
density of the inert gas, ncc; and on the mass of a single gas 
particle. The change of internal energy for an inert collision 
is simply the the change in kinetic energy after the collision: 



KE{T^, 



(14) 



where both KE{Tgas) and KE{Tg^y:a,in) are randomly drawn 
from Maxwell-Boltzmann distributions at the gas temper- 
ature Tgas and grain temperature Tgrain, respectively. We 
make the approximation that as the inert gas particle comes 
in contact with the cluster, it immediately attains the clus- 
ter's temperature, and rebounds with that dependency. 
Even though this is only an approximation, it ensures that 
over many collisions thermal equilibrium is established. 



where Cabs is the correction term resulting from the fact 
that clusters do not behave as blackbodies. This term is 
calculated by: 



a 



abs 



e(v,i)v dv 



V dv 



(16) 



where e[y, i) < 1 is the emissivity of the material, dependent 
on the size of the cluster, and the frequency of radiation 

(Draine & Lee 1984; Laor & Draine 1993). For small grains. 
Cabs severely reduces the rate of photon absorption as the 
cluster only significantly interacts with a few frequencies. 

The change in internal energy for a photon absorption 
is given by: 



A?7Abs = hv{T,^) 



(17) 



where v is randomly drawn from the modified blackbody 
distribution at the radiation field temperature Trad, given 
by: 



P(j/, i) oc 



(18) 



^rad — 1 



2.5 Photon Emission 

The rate at which the cluster emits photons is analogous to 
the rate of photon absorption: 



7i = 0.37(367r)i/^Ce^is ( —) ^££Z^i^/^ (19) 
except that it is evaluated at Tgrain and Cemis is given by: 

(20) 



f e{v,i)v'^ dv If dv 

Cemis = / t,^ / / fe^ 

J g^^grain" \ / J g ''^^grain \ 

The change in internal energy for a photon emission is: 

A?7Emis = -/w^(Tgrain) (21) 

where v is instead randomly drawn firom the modified black- 
body distribution at the grain temperature Tgrain: 



P{^,^) « 



g^'grain _ 1 



(22) 



2.6 Grain Temperature 

After each event, the internal energy of the cluster/grain is 
updated according to the above equations. Once the internal 
energy is known, the temperature of the cluster/grain can be 
calculated if the specific heat of the material is known. The 
new temperature is then fed back into the rate equations, a 
new event is computed and the cycle is repeated until either 
the cluster evaporates or grows into a stable grain. 



2.4 Photon Absorption 

The rate at which the cluster absorbs photons is derived 
from an integrated form of the Planck law, given by: 

= 0.37(36^)^/^Cabs f ^^^i^^ (15) 



3 MONTE CARLO CALCULATIONS 

The physical process described above can be implemented 
in a Monte Carlo computer code to follow the random os- 
cillations of the grain temperature and the grain fate under 
the various driving forces. Here and in the following, we ap- 
ply the theory described above to clusters of carbon atoms. 
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Figure 1. Temperature and size evolution of a grain initially 
composed of 1000 carbon atom immersed in a hydrogen-carbon 
gas with H/C = 100, Tgas = 3500 K, and with saturation 5 = 4. 



precursors of carbonaceous dust grains. We adopt the ab- 
sorption coefficients for graphite grains from Draine & Lee 
(1984) and Laor & Draine (1993), extrapolating their result 
to smaller grains when necessary. The other graphite prop- 
erties we adopt and relevant references are summarized in 
Table [21 

An example of the thermal history of a relatively big 
carbonaceous grain is shown in Figure (T) where a grain of ini- 
tially 1000 carbon atoms is left free to evolve in a hydrogen- 
carbon gas with number density ratio H/C = 100, temper- 
ature Tgas = 3500 K, and saturation 5 = 4. The grain is ob- 
served to grow as its temperature fluctuates with excursions 
as large as 200 K, or ~ 5 per cent of the grain temperature. 
Figure [2] shows instead the much more violent thermal his- 
tory of a small grain, injected in the same gas. In this case 
the temperature can fluctuate by more than 100 per cent. 
Even under such violent fluctuations, the grain eventually 
grows to a larger size. 



3.1 Grain Stability 

An important consequence of our analysis of the fluctuations 
of the grain temperature under the effect of quantized en- 
ergy exchanges with the gas and the radiation held is that 
the grain stability becomes a statistical concept instead of 
an exact concept. In the classical approach, when the Gibbs 
free energy is maximized, the grain is supposed to be in an 
unstable but exact equilibrium, able to remain at the sta- 
ble size forever. Any cluster smaller than the critical size is 
doomed to evaporate and any grain bigger than the criti- 
cal size will eventually grow. In our approach, two grains of 
the same size can be in different stability conditions (one 
stable and the other unstable) if they are at different tem- 
peratures. Moreover, a grain that is stable can become an 
unstable cluster without changing its size by absorbing a 
photon that raises its temperature and vice versa. A critical 
size that depends only on the gas temperature and satura- 
tion - reminiscent of the one of the classical theory - can 
be deflned as the size for which the probability of eventual 
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Figure 2. Same as Fig.[T]but for an initially smaller grain made 
by 10 carbon atoms. 



evaporation is the same as the probability of eventual growth 
into a stable grain. 

Here we discuss stability by finding the saturation at 
which a certain cluster becomes stable, according to the 
above definition. Unfortunately the photon absorption in- 
teractions become so numerous at low temperatures (Tgas < 
2500 K) that running the code until a statistically sound 
number of C atom attachments or detachments is observed 
becomes practically unfeasible. For that reason, in the fol- 
lowing we concentrate on the high temperature nucleation. 

Classically, the saturation at which a cluster with size i 
is stable, or critical, is given by: 



5- ^ gV 3 ; V PA- ; (23) 

Since Tgrain is not constant in our approach, we determine 
the stability saturation for a given cluster at a specific gas 
temperature by proceeding iteratively as follows. We start 
from a fairly large saturation (e.g. S = 10) and perform a 
minimum of 10 simulations in which the grain may evolve 
at that saturation. Each simulation is halted when the ini- 
tial cluster has changed its size by at least 10 per cen10. If 
all or most of the simulations resulted in grain growth, the 
estimated saturation was too large, and if all or most of the 
simulations resulted in grain evaporation, the estimated sat- 
uration was too small. We adjust our estimate of the critical 
cluster size and repeat the process, narrowing in on the satu- 
ration for which half of the clusters grow and half evaporate. 
In this way we can obtain a new stability curve, i.e., the lo- 
cation in the Tgas — S plane where a grain of a given size 
has the same probability of growing or evaporating. This is 
analogous to the critical cluster size in the classical nucle- 
ation theory. Thus, if the gas saturation is greater than the 
stable saturation the grain is more likely to grow. 

The stability plot for clusters with i = 1000 is shown 
in Figure [3] For all gas temperatures, the new stability 



^ For large clusters (j > 100) this is more than enough to de- 
termine the final fate of the cluster. However, for small grains 
we halt the simulations after a 70 per cent change of the cluster 
initial size. 
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Property Notes 

a = 1500 erg cm~^ Surface energy (Tabak et al. 1975) 

p = 2.23 g cm"'' Density of graphite 

mc = 1.9944~^^ g Mass of one carbon atom 

Eq. 3.3 of Guhathakurta & Draine (1989) Heat capacity 

Tbebyc = 420 K Debye temperature (Guhathakurta & Draine 1989) 

-844282 
13 ~T '■ — 

'n-x ea = — cm~^ Equilibrium gas density 

' "■-'grain 

^Bind.blk = 1.1831760"^^ erg Bulk binding energy of graphite 

Table 2. Graphite properties used in this paper. 
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Figure 3. Stability plot in the Tg^ain — 5 plane for a carbonaceous 
cluster with i = 1000 carbon atoms immersed in a hydrogen- 
carbon gas with H/C = 100. Both the solid line with square 
markers and the dash-dot line with star markers, show the re- 
sults of this paper. Actual calculations were performed only at 
the marked locations. Both lines show stability with temperature 
fluctuations, however, the solid line does not include an external 
radiation field, while the dash-dot line does, incorporating the 
heating effect of photons. The radiation field is assumed to be a 
black body at the temperature of the gas. For higher tempera- 
tures, radiation effects become negligible and the lines converge. 
The dashed line shows the stability for spherical clusters in the 
classical approximation in which the temperature of the cluster 
is constant and equal to the temperature of the gas. The thick 
solid line shows instead the stability line if the Einstein model is 
applied to the standard theory but temperature fluctuations are 
neglected. 

curve is lower than in the standard theory, generally al- 
lowing lower saturations in which clusters can grow. This 
is a result of three competing effects. At high temperatures 
(Tgas > 3000 K) the strongest effect is the Einstein model 
correction that makes the clusters much more resilient to 
evaporation. At lower temperatures, the grains are made 
more stable (even more than the Einstein correction predic- 
tion) by the fact that their temperature is lower than the 
gas temperature (Lazzati 2008). Both these stabilizing ef- 
fects are balanced by the effect of temperature ffuctuations 
that tend to make clusters unstable. The drop of the sta- 
bility line at low temperatures is expected to be reduced in 
the situation in which the radiation field is dense, so that 
interactions mediated by photons can also have a heating 
effect, as well as a cooling effect. In the case of a blackbody 
radiation field in thermal equilibrium with the gas, the drop 



Figure 4. Same as Figure|3]but for a smaller cluster with i = 10 
carbon atoms. Even with the effect of an external radiation field 
(stars) the Einstein correction dominates. 



is absent, as shown by the star symbols in Figure [S] This 
last condition may be relevant for dust production in a su- 
pernova explosions. 

Figure |4] shows instead the stability line for smaller clus- 
ters with j = 10 carbon atoms. In this case, the Einstein 
correction term has a big infiuence, allowing its effect to 
dominate over the entire temperature range that we investi- 
gated. Even the heating effect of a radiation field in thermal 
equilibrium with the gas does little to diminish the effect of 
the Einstein model correction (see the star symbols in Fig- 
ure |4l This steep drop makes determining stability for lower 
temperature more difficult. Also, smaller clusters are more 
susceptible to random fluctuations, inherent in our code, 
requiring more simulations, at least a minimum of 30, to 
determine the new stability curve. 

Another way to look at cluster stability and critical clus- 
ter size is to plot the attachment and detachment rates as 
a function of cluster size for given conditions of gas temper- 
ature and saturation. With the use of a Monte Carlo code 
like ours it is not straightforward to compute the average 
detachment rate for a given cluster size. That is because in 
a run, the cluster size varies (see Figure [1] and [2]) . In order 
to compute the average detachment rate from a cluster of a 
given size i we run the code artificially maintaining the clus- 
ter size constant. When an attachment happens, we compute 
the temperature that the i + 1 cluster would attain, but we 
keep the size equal to i. Analogously, when a detachment 
happens, we compute the temperature that the i — 1 clus- 
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Figure 5. Attachment and detachment rates from carbonaceous 
clusters for S = 2 saturation and gas temperature Tgas = 3000 K. 




T=3000K 



Einstein Corrected 
Temperature Fluctuations 
Classical Kinetic 
Classical Thermodynamic 



,Q-30 U d , , , , , 

^" 12345678 

Saturation Coefficient 

Figure 6. Nucleation rate as a function of saturation for car- 
bonaceous grains from a hydrogen-carbon gas with H/C = 100 
and temperature Tgas = 3000 K. 



3.2 Nucleation Rates 



ter would have but we maintain the size constant to i. This 
artificial forcing is equivalent to assuming that the temper- 
ature distribution of the size i cluster is not different from 
the temperature distribution of the i — 1 and i -I- 1 clusters. 

Figure [5] shows attachment and detachment rates from 
clusters containing between i = 3 and i = 10000 carbon 
atoms embedded in a saturation S = 2 hydrogen-carbon gas 
with H/C = 100 and temperature Tgas = 3000 K. The at- 
tachment rate is shown as a thick dashed line while detach- 
ment rates for different assumptions are overlaid. The clas- 
sical detachment rate is shown with a dash-dot line, yielding 
a critical cluster of ic = 350 carbon atoms. The dotted line 
shows the predicted rate if the Einstein model is considered 
but thermal fluctuations are neglected. The graph clearly 
shows that by itself the Einstein correction is too strong, 
making even the smallest clusters stable at a modest satu- 
ration. This would result in a catastrophic nucleation of very 
small grains. The detachment rate with thermal fluctuations 
included is shown with a black line. Interestingly, the at- 
tachment and detachment rates almost trace each other for 
the small, unstable clusters. This results from the fact that 
the biggest thermal fluctuations are due to the attachment- 
detachment process, since the binding energy is typically 
much larger than the photon and collision energies. As a 
consequence, an attachment results in a sharp increase in 
the cluster temperature that increases the detachment rate 
by a large factor, increasing the probability of a detachment, 
which in turn cools the cluster significantly. This locking of 
the attachment and detachment events is particularly evi- 
dent for pure gases and for fcs = 1, in which every collision 
of the cluster is with a monomer that attaches to the clus- 
ter. Less pure gases show this tracking to a lesser extent 
(typically only for small clusters). In terms of the critical 
cluster size, we see that the critical cluster is smaller than 
the one predicted classically, but not by a huge factor. The 
complete theory with fluctuations predicts a critical cluster 
of approximately ic ~ 100 carbon atoms. 



Once the average rates out for all sizes between i = 3 and 
approximately twice the critical size have been computed, 
the nucleation rate can be calculated with a modiflcation of 

Eq.m 



J = 



nx 



2ic 



J=3 



(24) 



Since the detachment rate is always signiflcantly less than 
the attachment rate for i = 2, due either to linearity of the 
molecular bond (thus a dramatic change in the internal en- 
ergy) or from a drastic drop associated with the Einstein 
correction term, we realize that we cannot clearly define a 
nucleation rate directly from the gas phase. Instead, we as- 
sume that all of the gas particles have formed i = 2 clusters, 
and compute J using half of the number density of the gas. 
The extreme robustness of the i — 2 clusters is probably due 
to the failure of the capillary approximation for very small 
clusters, as discussed below. 

Figure [B] shows the rates computed with our Monte 
Carlo code compared with rates obtained by other methods. 
The dash-dot line shows the rate obtained from the thermo- 
dynamic approach (e.g. Kozasa et al. 1989, 1991; Todini & 
Ferrara 2001; Nozawa et al. 2003, 2008, 2010), while the 
dashed line shows the rate obtained from the kinetic theory 
under the assumption of thermal equilibrium between the 
cluster and the gas. A solid line shows the nucleation rate 
if the Einstein model correction is taken into account, but 
thermal fiuctuations are neglected. As expected, this rate 
is very high since the correction makes clusters of all sizes 
stable. Finally, a solid line with circles shows the result with 
temperature fluctuations. Since, as we saw before, the com- 
bined effect is that of making the clusters more stable, the 
nucleation rate is significantly larger than the classical re- 
sult, especially at moderate saturation. Figure [7] shows the 
same results for a hotter gas with Tgas ~ 4000 K. Finally, 
Figure [5] shows a contour plot of the nucleation rate as a 
function of both gas temperature and saturation. It appears 
that the nucleation rate is highly dependent on temperature 
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Figure 7. Same as Figure [6] but for a gas at temperature Tgas ■ 
4000 K. 



Figure 8. The log of the nucleation rate for carbonaceous grains 
from a hydrogen-carbon gas with H/C = 100. 



at high saturations, but becomes increasingly dependent on 
saturation as the threshold 5* = 2 approaches. 

A final, important result, if not the most important 
result of this work, is shown in Figures |9] and 1101 The 
figures shows the nucleation rate of carbonaceous grains 
as a function of the carbon concentration in the gas, for 
10 < H/C < 300 and for various saturation values. Ac- 
cording to the classical theory, both thermodynamic and ki- 
netic, the nucleation rate should only depend on the partial 
pressure (or partial density) of the compound that partici- 
pates in the nucleation, carbon in our case. Instead, we see 
that if thermal fluctuations are relevant, the nucleation rate 
depends quite strongly on the concentration. The depen- 
dence of nucleation rate on concentration is complex and 
not straightforward to interpret. For low saturation, the nu- 
cleation rate decreases as the relative density of the carrier 
gas increases, i.e., a pure carbon gas nucleates much faster 
than a mixed gas (Figure [U]). At intermediate saturations, 
the fastest nucleation is observed at intermediate concen- 
trations {H/C ~ 50 for 5* = 4, lower panel of Figure [10)) . 
Finally at high saturations, the nucleation rate increases at 
low concentrations (or high H/C values). The trend seen 
in Figure |9] at low saturations is easy to explain. The pres- 
ence of carrier gas particles creates a thermal contact be- 
tween the cluster and the gas that established thermal equi- 
librium without requiring growth or evaporation, breaking 
the strict growth^heating^evaporation— ^-cooling— ^-growth 
behavior outlined above and making the clusters more vul- 
nerable. The trends observed in Figure [TO] are instead more 
difficult to fully understand. Inspection of the rate plots 
(analogous to Figure [5| suggests that the increase of the 
nucleation rate is due to the behavior of the grains just 
above the critical size rather than to the attachment and 
detachment rates from sub-critical clusters. 



4 DISCUSSION 

We have presented a model for the physics of the nucleation 
rate of dust grains from the gaseous phase that includes 
the role of thermal fluctuations of the clusters as well as 



corrections for the detachment rates from nanoscale clus- 
ters. Our implementation is based on the kinetic theory of 
nucleation (Becker & Doring 1935, Kashchiev 2000), but in- 
corporates some fundamental changes in the way in which 
the thermal balance of the forming clusters is considered. 
First, we included a modification to the detachment rates 
from small grains, following the Einstein model developed 
by Guhathakurta & Draine (1989). This is a stabilization 
effect that takes into account the fact that very small grains 
at low temperature need to use a large fraction of their total 
internal energy to brake a bond and eject a monomer. As a 
consequence, the ejection rate does not scale linearly with 
the surface but becomes very small for small grains. In addi- 
tion, we consider the fiuctuations of the grain temperature 
as the grain exchanges energy with the gas and the radiation 
field. Since all the energy exchanges are quantized (through 
photons or collisions), the temperature of the clusters does 
not attain an equilibrium value but rather continuously fluc- 
tuates around it. We find that the combination of these two 
effects results in an increased tendency of the clusters to 
grow, and as a consequence, yields higher nucleation rates. 
We also find that the nucleation rate depends on the con- 
centration of the nucleating compound in the inert carrier 
gas, with higher rates for higher concentrations, especially at 
moderate saturation. This is a new result that has not been 
predicted by any other theory and can potentially become 
a crucial test for this model. 

Our study has, like most nucleation studies, some im- 
portant limitations that deserve some discussion and could 
sizably affect the numerical values of the critical clusters and 
nucleation rates. 

• First, we were able to numerically implement the theory 
only for high temperatures, for which the rate of particle in- 
teractions (cluster-monomer or cluster-carrier gas) outnum- 
ber the radiation interactions (carrier-photon). We cannot 
therefore give any result for the interesting cold tempera- 
ture regime, where supercooling of the cluster can happen, 
as discussed in Lazzati (2008). A possible solution for the 
low temperature regime is to treat radiation interactions 
as continuous and concentrate on the quantized nature of 
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Figure 9. Nucleation rate as a function of the carbon concen- 
tration in the gas phase for a iiydrogcn-carbon gas at saturations 
1.2 and 1.5 and temperature 4000 K. For low saturations, higher 
H/C ratios drastically lower nucleation rates. 

the particle interactions only. Such implementation is under 

study. 

• More fundamentally, we have been implicitly adopting 
the capillary approximation, i.e., we have assumed that the 
surface energy of graphite clusters is the same of the one of 
big graphite chunks and we have assumed that the concept 
itself of surface energy applies to clusters as small as a hand- 
ful of carbon atoms, for which a surface is a hard concept to 
define. It is likely that this results in an overestimate of the 
binding energies of small clusters, ultimately resulting in an 
overestimate of the nucleation rate. 

• Another potential overestimating factor is that wc as- 
sumed the sticking coefficient to be fcs = 1, i.e., all incoming 
monomers that collide with the cluster stick to it and make 
the cluster grow. The nucleation rate is directly proportional 
to the sticking coefficient and, since ks < 1, our nucleation 
rates are overestimated if the real coefficient is different from 
1. In addition, reality is likely more complex, ks is probably 
dependent on the number of atoms in the cluster. Clusters 
that fill all the bonds of carbon atoms (such as fuUerenes, 
for example), arc likely to have very small values of ks since 
it is hard to make tlieiii grow beyond their closed structure. 
An i-dependent value of ks would, again, result in a decrease 
of the nucleation rate. 

• We have neglected the possibility of cluster-cluster in- 
teractions. Even though these should not be numerous, they 
can be relevant for a particularly pure gas of the nucleating 
compound. 

• We have assumed that even the smallest clusters emit 

continuum radiation. It is instead likely that they will emit 
radiation in bands becoming more and more narrow as the 
size decreases. 

• We have assumed that all the energy released as a bond- 
ing of a new monomer happens is transformed into internal 
energy of the cluster/grain. However, it is possible that some 
of that energy is carried away by a photon. If this were true, 
smaller temperature fluctuations would be seen in conse- 
quence of attachments of new monomers. 

Despite all these limitations, we believe the two main 
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Figure 10. Nucleation rate as a function of the carbon concen- 
tration in the gas phase for a hydrogen-carbon gas at saturations 
4.0 and 8.0 for 4000 K. There appears to be a transition such that 
higher saturations instead allow higher H/C ratios to increase nu- 
cleation rates. 

conclusions of this work hold. First, thermal fluctuations are 
important in the nucleation process. Second, if this is true, 
the nucleation rate should depend on the concentration of 
the nucleating compound. The difference of the nucleation 
rates resulting from different theories should warn us that a 
lot has still to be learned in this fascinating fleld and conclu- 
sions based on the blind application of classical nucleation 
rates should be taken with care. 
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